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Abstract. We resolve several longstanding problems concerning the stability 
and the absence of multi-particle binding for N > 2 polarons. Frohlich's 
1937 polaron model describes non-relativistic particles interacting with a scalar 
quantized field with coupling y^, and with each other by Coulomb repulsion of 
strength U. We prove the following: (i) While there is a known thermodynamic 
instability for U < 2a, stability of matter does hold for U > 2a, that is, the 
ground state energy per particle has a finite limit as N —> oo. (ii) There is no 
binding of any kind if U exceeds a critical value that depends on a but not on 
N. The same results are shown to hold for the Pekar-Tomasevich model. 



Frohlich's large polaron jl] is a model for the motion of an electron in a polar 
crystal and it is also relevant as a simple model of non-relativistic quantum field 
theory. Consequently there is a huge literature, both experimental and theoretical, 
devoted to its study. See, e.g., [12l[32l|3Tl|27l|T] and references therein. Our concern 
here is with the binding or non-binding of several polarons: whether the ordinary 
Coulomb repulsion among the electrons can, if strong enough, prevent the binding 
that would otherwise be created by the electric field of the polar crystal. We are 
also interested in the stability of matter, i.e., whether the energy of N polarons is 
bounded below by a constant times N even when there is binding. 

In this model the single polaron, which is one non-relativistic electron interacting 
with a phonon field, has the Hamiltonian 



This Hamiltonian acts in the Hilbert space L^(R'^) (g) where is the bosonic 
Fock space for the longitudinal optical modes of the crystal, with scalar creation 
and annihilation operators a^(fc) and a(fc) satisfying [a(fc), a^(fc')] = S{k — k'). The 
electron momentum is p = — iV, the phonon field energy is 



1. Introduction and main results 



(1.1) 




(1.2) 



and the interaction of the crystal modes with the electron is 



0(x) = ^ / ^ (e'^-a(fc) + e-^-at(fc)) 
V27r J«3 \k\ 



(1.3) 



Date: April 27, 2010. 
© 2010 by the authors, 
purposes. 



This paper may be reproduced, in its entirety, for non-commercial 



2 



R.L. FRANK, E.H. LIEB, R. SEIRINGER, AND L.E. THOMAS 



with coupling constant a > 0. (In another frequently used convention a is replaced 
by a/y/2.) The ground state energy E'-^^{a) is the infinium of the spectrum of H'-^\ 
Because of translation invariance, E^^\a) cannot be expected to be an eigenvalue, 
and indeed it is not; this was proved in [11] using methods developed in (lOj . 

A noteworthy feature of the phonon field energy Hf is its flat dispersion relation, 
i.e., there is no non-constant function a;(A;) in the integrand of (jl.2p . The energy 
of infrared phonons does not go to zero as |fc| — >■ 0, while the ultraviolet energy 
is finite when |fc| — )• oo. If we tried to minimize the energy in a nai've way by 
completing the square, we would end up with a Coulomb- like J dk\k\~'^ self-energy 
of the polaron. The non-integrability for large k would lead to a divergent self- 
energy, but this divergence is actually mitigated by the electron kinetic energy 
and the uncertainty principle. While the single polaron has finite energy, another 
problem remains for the many-polaron system; the energy is finite, but stability 
of matter will not hold unless a sufficiently strong Coulomb repulsion among the 
electrons is included in the Hamiltonian. 

The Hamiltonian for N electrons is 

N 

= E (P' - \^<^(^0) +Hf + U Vc{X) (1.4) 
1=1 

with X = {xi, . . . ,xn) e and 

We impose no symmetry restrictions on the electrons, which means that our lower 
bounds apply equally to bosons or fermions or particles with no symmetry restric- 
tions (boltzons). The Hilbert space is then i^(R'^^) ^T. Particle spin is irrelevant 
for our results and is ignored. Physically, the parameter U is the square of the 
electron charge, and it satisfies C/ > 2a Nevertheless, we will consider all values 
U>0. 

The ground state energy of Hjj''^ is denoted by E^^\a), and the binding energy 

is defined to be AElj^\a) — NE'^^\a) — E\f\a). We will prove three theorems 
about these quantities. They were previously summarized in an announcement [5]. 
Our main goal is to find conditions on U and a such that no binding occurs, i.e., 
A.E\f\a) = 0. Of particular physical interest is the case N = 2 (bipolaron). 

Theorem 1 (Absence of binding for N polarons). For given a > there is a 
finite Uc{a) > 2a such that 

AElf\a)^0 forallN>2 (1.6) 

whenever U > Uc{a). 

Our proof is constructive and gives an explicit upper bound on Uc{a); see the 
discussion at the end of Section H) This bound on Uc{a) is linear in a for large a, 
which is the correct behavior. Presumably, the true Uc{a) behaves linearly even 
for small a, but this remains an open problem. 

For the bipolaron, N = 2, our proof is simpler and yields the sharper result that 
the critical Uc{a) indeed obeys a linear law: 
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Theorem 2 (Absence of binding for bipolarons). Let N = 2. For some 
constant C < 26.6, 

AEl^\a)^0 (1.7) 

whenever U > 2Ca. 

The optimal constant C in Ttieorem [2] is presumably much closer to 1 than the 
bound we derive. It is not equal to 1, however. In the strong coupling limit a — >■ oo, 
binding occurs if f7 < 2.3a [331 IS]- For small a, on the other hand, variational 
calculations [331 H] suggest that bipolaron binding does not occur at all for any 
U > 2a. The proof of this remains an open problem. 

On the other hand, if binding does occur, we would like to know how the binding 
energy depends on N; in particular, is there stability of matter, in the sense that 
AElf \a) < C(C/, a)N for all N7 

This linear bound, if it exists, implies the existence of the thermodynamic limit 

lim N-^E^r^\a). 

The proof is a simple consequence of the sub-additivity of the energy, i.e., 

E[r^''\a)<Elf\a) + E^f\a), (1.8) 

which follows from the fact that one can construct variational functions in which 
N electrons are localized on the earth and M behind the moon [13] [20l Sec. 14.2]. 

The proof of the linear lower bound is far from obvious; indeed, it is not always 
true! Griesemer and M0ller [13j recently proved that when U < 2a, there are posi- 
tive constants (depending on U and a) such that —ciN'^/^ > E\f\a) > —C2N'^^'^. 
This result holds for particles, like electrons, that satisfy Fermi statistics. If the 
electrons were bosons, the result would be even worse, —N^, as a similar analy- 
sis shows. This is the same behavior as that of gravitating particles in stars [20l 
Ch. 13]. In the opposite regime, C/ > 2a, [13] shows that E[f \a) > ''C{U,a)N'^, 
independently of statistics, where C{U,a) — oo as C/ \ 2a. (In their convention 
the dividing line is U = V2a.) 

Judging from the physics of the model, it is reasonable to suppose that there is 
a linear law as soon as ?7 > 2a. This we are able to prove. 

Theorem 3 (Stability for U > 2a). For given U > 2a > 0, N~'^E^\a) is 
bounded independently of N . 

Our lower bound on N~^E^\a) goes to — oo as 17 \ 2a, but we are not 

claiming that this reflects the true state of affairs. Whether limAr_j.oo N~^E'2^\ct) 
is finite or not remains an open problem; see, however, the discussion of the strong 
coupling limit below. 

For U in the range 2a < U < Uc{a), there are bound states of an undetermined 
nature. Does the system become a gas of bipolarons or does it coalesce into a true 
iV-particle bound state? If the latter, does this state exhibit a periodic structure, 
thereby forming a super-crystal on top of the underlying lattice of atoms? This is 
perhaps the physically most interesting open problem. 

1.1. The strong coupling limit. There is a non- linear differential- integral varia- 
tional principle associated with the polaron problem, which gives the exact ground 
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state energy in the limit a — >■ oo. This variational problem was investigated in de- 
tail by Pekar [28] . Pekar and Tomasevich (PT) [29] generalized it to the bipolaron, 
and the extension to A^-polarons obviously follows from [55^. 

The PT functional is the result of a variational calculation and therefore gives an 
upper bound to the ground state energy E^\a). In order to compute Hjj'^'^), 
one takes a 5' of the form ^ (g) $ where 'ip G i^(R^^), $ G J", and both and $ 
are normalized. For a given -0 it is easy to compute the optimum $, and one ends 
up with the functional 



■dxdy, (1.9) 



i<j 



where dX = IlfeLi dxk, and 



N 



p^(x) = > / \iP{xi,...,x,...,xn)\ dxi- ■ ■ dxi- ■ ■ dxM (1-10) 

~[ JR3(«-1) 



with X at the i-th position, and dxi meaning that dxi has to be omitted in the 
product J^^2 dx].. The ground state energy is 

S\f\a)=ini\v\f\^] : / m-'dX = l\ . (1.11) 



■■u V"; — ^^ii ^ I u 

Hence the variational argument above gives the upper bound 

i?(f)(a)<£r(a)= 4)1(1) (1.12) 

(The equality follows by scaling.) For = 1 this upper bound is due to Pekar; nu- 
merically, one has f ^""^^(a) ~ — (0.109)q;^ Moreover, the minimization problem 
for E^^\a) has a unique minimizer (up to translations), see [IB) . 

The upper bound for iV = 1 was widely understood to be asymptotically exact 
for large a. A proof of this was finally achieved by Donsker and Varadhan [Bl, using 
large deviation theory applied to the functional integral discussed below. Later, this 
fact was rederived in [23] by operator methods, and it was shown that the error 
was no worse than a^/^ for large a. 

The fact that for fixed ratio v — U / a >Q 

lim a-^E\^\a) = 8i^\a = 1) (1.13) 

for N — 2 and any i' > was first noted in ^S]. This is also valid for arbitrary N. 

It follows from the limiting relation (|1.13p , together with the fact that our bound 
on Uc{a) is linear for large a, that our three theorems about E^\a) transfer to 
the same theorems about E\f \a) which we state next. 

Corollary 1 (Stability and absence of binding for the PT functional). 

1. Stability holds for U > 2a, that is, for any v > 2 there is a constant C(j^) such 
that 

£u^\a) > -C{v)a^N for all N > 2 (1.14) 
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whenever U — va > 0. 

2. There is no binding if U/a is large enough, that is, there is a finite Vc > 2 such 
that 

£lf\a)=N£^^\a) for all N > 2 (1.15) 

whenever U > t'cO;. 

We note in passing that the critical U in the PT model depends linearly on a. 
This follows by scaling. 

In Section [5] we will give a direct and easier proof of Corollary [1] that has no need 
of the functional integral machinery and leads to better constants. In particular, 
with our new proof we find that 

r(0.0280).3/(._2) if2<.<3, 
~ [0.755 if 1/ > 3. 

For fermions, a stronger result than Corollary [T] was obtained in |13] . They prove 
that at critical coupling U ~ 2a one has 

£^^\a,q) > -{0A61)^a^q^/^N for all iV > 2 and ah 1 < q < , (1.17) 

where sjj'^ (a, q) is the infimum of vlf"^ [^p] restricted to fermions with q spin states 
{q = 2 for electrons). Recall that the single polaron energy is £^^\a) ~ — (0.109)a^ 
[25] . For completeness we repeat the short proof of (|1.17|) in Section [Sj 

1.2. Previous results on the polaron ground state energy. In addition to the 
large a asymptotics just mentioned, there have been several other rigorous results 
on the ground state energy of the (single) polaron, some of which we will use here. 

(i) One of the earliest results was the variational calculation of E'^^^(a) for 
small a by Gurari and by Lee, Low, and Pines [HI [TCI [T7] which leads to 

^;(^)(a)<-a foraUa. (1.18) 

(ii) A lower bound, which validates the conclusion that E'^^^ (a) ^ —a as a — > 0, 
was obtained in [l^. They prove that 

E'^^\a)>-a-\o? for ah a. (1.19) 

(To derive pTTOl use p = 1 + 2q/3 in [24, Eq. (24)].) We note the fact that 
()1.19|) has the correct power law behavior for both small and large a. 

(iii) The functional integral formulation: The large time behavior of the heat 
kernel exp{-TH^^^) gives us the ground state energy as 

E[j^\a) = - ^lini In {Jexp ( - Ti/^f ^)) , (1.20) 

where (•) denotes the expectation in a suitable state, i.e., a normalized 
vector in the Hilbert space. Since the phonon operators can be realized 
as the coordinates of a quantum-mechanical harmonic oscillator (one for 
each value of fc), we can apply the Feynman-Kac formula for the evaluation 
of (exp(— Tff^^'')). Since the harmonic oscillator coordinates appear only 
linearly and quadratically in the exponent, they can then be integrated out 
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explicitly. One obtains a formula due to Feynman [7], see also [30l Sec. 5.3] 
for a careful discussion. The conclusion is that 

Elf\a) = lim lim T-MnZgXT) , (1.21) 



where 



ZuHiT)-^ I dxi--- I dxN I dW:i'^{uji)---dW'^'^{LJN)xBniu}i)---XBn{uJN) 

I f rfse-i^i ^ dt 

X exp a / 2^ 



/ dxi ■ ■ 


j dxN 


Ibr 


J Br J 



\Jo \uJ^{t)-UJj{t + s)\ 

xexpf-[/V/ - — — TTT 1 ■ (1-22) 



i<j 

Here dW^ denotes the Wiener measure of closed Brownian paths in 
with period T starting and ending at x. Moreover, Bji denotes the ball 
centered at the origin of radius R, and the characteristic function of the 
path XBri('^j) is 1 if '^j stays inside the ball Br for all times, and zero 
otherwise. The argument of iOj{t + s) is understood modulo T. 

2. The two-polaron problem: Absence of Binding 

We first consider the special case N = 2 and prove Theorem [2] It is convenient 
to structure the proof in three steps. 

Step 1. Partition of the interparticle distance. We choose a quadratic partition 
of unity (IMS localization [H Thm 3.2]) and localize the particles according to their 
relative distance. (In the A'^-particle case later on, we will localize with respect 
to the nearest neighbor distance, which, for = 2, is the same as the relative 
distance.) This kind of localization is one of the principal novel features of our 
analysis. 

In order to construct this partition, we pick some parameters 6 > 1 and £ > 0, 
and let 

krt<£/b, 



cosf^ for^<t <6^, ' ' 

fort>be. 

For j > 1, let ipj{t) ip{b'^-H), and for j = 0, let 

r 1 for t < e/b, 

Mt)-^{ cosf^ torl/b<t<e, (2.2) 
[ for t > £ . 

Then 

^j{tf = 1 for alH > . (2.3) 
Using the IMS localization formula, we can write, for any wave function -0, 



iorl/b<t<i, 



X2\ 



V',)=:E^^-|I^^H' (2-4) 
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with ipj{xi,X2) — '4>{xi,X2)'^j{\xi — X2\) and with numbers Cj (depending on ijjj). 
Our goal is to prove that Cj > 2E'^^\a) for all j if U > 2Ca. If this is indeed the 
case, then the right side of exceeds 2E^'^'>{a) J^j llV'jiP = '2E'^^'> {a)\\^p\\^ , which 
is the assertion of the theorem. 

For our bounds we shall use the fact that on the support of ifjQxi — X2\), the 
localization error is dominated by 

^2 



k>0 



\^'k{\xi -X2\)\ 



< 



Mi-i/b)^ 



1 ifj = 0, 

fe2(l-j) if j > 1 . 



(2.5) 



Moreover, on these supports, we shall bound the Coulomb repulsion from below by 

U 



Xi 



X2\ ~ I 



for all J > . 



(2.6) 



It is clear from (|2.5p and p.6p that by choosing U large enough, we can dominate the 
negative localization error by a part of the positive Coulomb term. What remains 
is to dominate the polaronic attraction by the remainder of the Coulomb repulsion. 
For this, we distinguish between the cases j > 1 and j — 0. 

Step 2. The case j > 1; Energy estimate for separated particles. We further 
localize each of the two particles to its own ball of radius VL for some parameter 
i > 0. This will entail an additional localization error. Concretely, let 

1 ( sin(7r|a;|) 




/2n\x\ 



for |a;| < 1 ^ 
for bl > 1 . 



(2.7) 



With 



and note that J dxx{xY — 1 and J (ix|Vx(x)p — 

i^],ux,u2{xi,X2) = i}j{xi,X2){y ly'^xib^' {xi - ui)/L)x{h^^{x2 - U2)/L) (2.. 
we have, by a continuous version of the IMS localization formula. 



^3 



4^)-2^|^',(l^i-^2| 



fe>0 



= / dui I du2 ( tpj,ui,U2 



(2) 



> 



du2 ( 1pj,Ui,U2 



-2^l'Pfe(ki-a;2| 

fc>0 

^ \2{i-i/b)^ 



2iivxir 



2^ 
L2 



(2.9) 



The latter inequality comes from (|2.5p . Note that since |a;i — X2I > V^^l on the 
support of ipj , the wave function V'i.ui ,U2 is non-zero only if the two balls of radius 
V L centered at Ui and U2, respectively, are separated at least a distance 

d>V-^£~AVL. (2.10) 

We choose this to be positive by requiring that L < £/{Ab'^). 

Lemma 1. Assume that ■0 is normalized and supported in Bi x B2 where Bi and 
B2 are disjoint balls of some radius R, separated a distance d. Then 



(V-liyf 1^) >2£;W(a)-^ 



2a 



(2.11) 
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This lemma will be proved in Subsection 12.11 It is an easy consequence of the 
functional integral representation of the ground state energy. 

We apply inequality (P?TT]) to ((^ . Using the bounds (PTTU)) on d and on 
the Coulomb potential, we conclude that Cj is bounded from below as 

This last expression is > iE'^^^a) for all j > 1 if and only if 

(the j > 1 condition) . (2.13) 



2al Mtt^ 2tt^ 

U > 



e/b^-AL 2{e-e/b)^ m 



Step 3. The case j — 0; Energy estimate for neighboring particles. Because of 
751) and (|2.6p we have the lower bound 

eo>i^r(«)-2(I^ + 7' ^'-''^ 

where Eq (a) denotes the two-polaron energy in the absence of Coulomb repulsion, 
i.e., for U ~ 0. The following lemma compares this energy with 2i;(i)(a). 

Lemma 2. For all a > 0, 

(a)>2i?W(a)-Ia2. (2.15) 

Also this lemma uses the path integral formulation and we defer the proof to 
Subsection 12.11 At this point we will utilize it to conclude the proof of Theorem 
[21 The constant 7/3 in (|2.15p is certainly not optimal, and an improvement would 
lead to a better constant in Theorem |2l 

It follows from ([2T4)) and ((21^ that eq > 2E'^'^\a) if 



U > W 



3 2{l-llbf 



(the j = condition) . (2-16) 



Numerical evaluation shows that the two conditions (|2.13p and (|2.16p on U can be 
satisfied for an appropriate choice of 6, £, and L if [/ > 61 a. (Choose b — 1.2, 
i = 22.8 and L = 0.142 ^.) For U satisfying these conditions, each ej > 
2E'-^'>{a). This completes the proof of Theorem |2l with the bound on the constant 
C < 30.5. □ 

In order to improve the bound on the constant C of Theorem [21 we replace 
Lemma [H by the following alternative bound. 

Lemma 3. Under the same assumptions as in Lemma{^ 



|V')>2£;W(a)-^V 



2a \ 16aR 
\xi-X2\ I ~ 7r2d(rf + 4i?) ■ ^' ' 

In the appendix we provide a proof of this lemma using the Rayleigh-Ritz vari- 
ational principle. The proof is certainly not easier than the one of Lemma [1] using 
the functional integral method, but it may be of use in other applications where a 
functional integral approach is not as convenient (or available). The method does 
point up the utility of localizing the phonon field about the respective particles - 
in this case, to half-spaces each containing a particle. 
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We apply the bound (f2l7l) to ([2?9l) . with R = WL and d satisfying (|2?T0l) . and 

conclude that 

(2.18) 

This expression is > 2E^^'> (a) for all j > 1 if and only if 



16aLb^ bin^ 2TrH 

- " TT^e/b^ - 4L) ^ 2{i - ^/6)2 + 



(2.19) 



Eqs. ((2l6| and ((2l9| are satisfied for U > 53.2 a with the choice b = 1.23, £ = 
19.7a-i and L = 0.15£. 

2.1. Some uses of the path integral. Lemmas [1] and [H used in the previous 
subsection, will be proved here. 

Proof of Lemma [H We use a Feynman-Kac representation similar to (jl.2ip . It 
implies that the infimum of the left side of (12. lip over all tp with the required 
support properties equals 

- lim ilnZB„B,(T) (2.20) 

T—^oo 1 



where 



Zb„bAT)-= I dx, I dx2 I dW^^{uJi)dWj^{u:2)XBA<^l)XBA^2) 

{ r dse-l'^l /"^ dt \ 

^'"P rl^^/^ LTTTY^TTTrr^ I • ^^.21) 



'^Jo \^^^it)-UJj{t + s) 



Here dWj^ denotes the Wiener measure of closed Brownian paths in M.^ with period 
T starting and ending at Xj, and XBj i^j) is 1 if i-^j stays inside the ball Bj for all 
times, and zero otherwise. Since |wi(i) — W2(^ + s)| ^ d for all t and s we see that 
Zbi.B2{T) is bounded from above by 

Replacing XBji'^j) by its upper bound 1, we deduce inequality (12. lip . □ 

Proof of Lemma [3 Application of the Cauchy-Schwarz inequality in the path inte- 
gral (IL^ yields 



Z'oliTr< // dW^iu;,)dW^{u;2) 



X ex (2a f V r rf^ \ 

dW^iuJl)dW^{i02) 

X exp 4q / / - — — ; , 2.22 

A 2 Jo M) - L02it + .s)\ ^ ' 
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where / diyj(cj) is short for Jg^ dx J dWj(w)xSj,, (t^)- The first factor on the right 
side equals the square of 

dW^iC(-)exp (2a f r ^\ ] , (2.23) 

\ Jw 2 Jo \uj{t) - uj{t + s)\ J 

which, in turn, is the one-polaron expression with a replaced by 2a. Using Jensen's 
inequality, we bound the second factor from above by 

Since closed Brownian paths are invariant under time reparametrization, the latter 
integral does not actually depend on s, and hence (I2.24p equals 

This functional integral represents two particles in the ball Bji interacting via an 
attractive Coulomb potential —'ia/\xi — X2\- This is like the positronium Hamil- 
tonian whose ground state energy equals —2a^ in the limit i? — > cx). Summarizing, 
after taking the T — >■ 00 limit we find that 

4^^(a) >£;(i)(2a)-a2. (2.26) 

To finish the proof, we use the bounds (|1.18p and (|1.19l) . which imply that 

E'^^\2a)>2E^^\a)~\a^ . (2.27) 
This, together with (H^, proves ((^1^ . □ 

3. The A^-polaron problem: Thermodynamic Stability 

We now consider the case of general N and prove Theorem [3] We start by 
localizing particles in balls in order to reduce the problem to a local one. We use 
the sliding technique introduced in (see also [51^ ) . Pick an even and real- valued 
function % with compact support, normalized by J = 1, and w > large enough 
such that the function 



/(^) = ^(l-e-l-lx*x(^)) (3.1) 



is positive definite. (The symbol * means convolution.) The existence of such an 
Lo for smooth enough x was shown in [3", Lemma 2.1]. For any operator- valued 
function p(x), 

dx dy [Y, S{x - X,) - pixA fix - v) ( ^ 5{y - x,) - p{y) \ > . (3.2) 



We apply this to 



p(x) = I dk \k\e'''''a(k) (3.3) 

' (2^)2^2^,' I I w \ J 



and obtain the bound 

V -. --^/7^y^(t>ix,) + Hf>~aNLu + 2a dzl^{z). (3.4) 

l<^<j<JV l^'^^^l .=1 
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Here, 



l<i<j<N ' * ^' 



1 r f fy~^\^i~y\ 

-2Yxz{x.i) J dy _ xz{y) {p{y) + p{yy) 



i=l 



1 ff g-'^k-yl 
+ 2JJ d^dyXz{x)p{xy |^_^| P{v)Xz{y) (3.5) 

where we denote Xz{x) = x{^ ^ •2;). We also note that 

= J dzpxlp^ J dzXzP^Xz- j da;|Vx(a;)p, (3.6) 

and thus 

Hlf^ > j dzHz + {U ~ 2a) Vc - aNuj - ^ J |Vx(a;)P (3.7) 



where 

N 

I 



1 ^ 

2Y{Xz(xt)pfxzixi) +PiXz{xi)'^Pt) +2alu,{z). (3.8) 
1=1 

The Hamihonian Hz is concerned only with the particles in the support of Xz', 
similarly for the phonon field, p{y) enters only for y in this support. Moreover, 
Hz commutes with Uz — X]j=i ^^(^^Oi the number of particles in the support of 
Xz , where 9z denotes the characteristic function of the support of Xz ■ We can thus 
look for a lower bound on Hz in a fixed sector of Uz particles. We will prove the 
following lower bound. 



Lemma 4. With [t]+ = max{t, 0}, 



N 



2 



ff.>-a[4an.-^]+5^X.(x.)'-^ (^^^llxlloo + llVxIbj • (3.9) 

Proof. Pick some A > a;. Applying p.2p with the positive definite function f{x) = 
|a;|~-^(e~"l^l — e""^'"^'), we have 

1 ^ 

U^) > Ia{z) - - {A - u;)Y,Xzix.)' . (3.10) 

i=l 

The last term represent the 'self-energy' terms. We will choose A proportional to 
Uz, hence this term is of the order 71^. For the remaining terms, we complete the 
square and write 

dy-: — -^—rXziy) {p{y) + p{y)^) 
\Xi y\ 

= {Xz{xi)pi - Az{xi)) (ptXzixi) - Azixi)'') 
+ (PiXzixi) + Az{xi)^) {xz{xi)pi + Az{xi)) 

- 2Az{x,)^ Az{x,) - [Az{x,),Az{x,)^] , (3.11) 
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where Az{x) is a vector operator with three components, 

:= -J y dyx.{y)p{y) J rffc ^2(fc2 + A2) ■ (3-12) 
The Schwarz inequahty apphed to this last equation shows that 

,-A|to-a| r, 2 



f ^-J^\w-y\ f 1 

Az{xyAz{x)< J dwdyxz{w)p{wy |^ _ ^| P{y)Xz{y) J dk^-^^ 



A2) ■ 
(3.13) 



Moreover, 



3 

^ E (a"'/'IIxIIooI|v/,||2 + llVxIbll/.lloo)' (3.14) 

with 

a:, l-(l + |x|)e-l"l 

/i(a;) = -Tj ri^ 7 x = (xi,a;2,X3) . (3.15) 



We choose A such that 



^fe fc2(fc2^A2) -^<"^ (3-16) 



TT 



which then assures that the last two terms of equation p. lip are relatively bounded 
by the field energy terms, i.e., the last term of Iu:{z) in p. 51) . More precisely, 
with the choice A = max{a;,4anz} Lemma 4 follows, using the facts ||/i||oc = 1, 
II ^/i II 2 = 27r/3, and the commutator bound p. 141) . this bound contributing the 
last term on the right side of p.9p . □ 



We now complete the proof of Theorem 3. If we insert bound p.9p into p.7p . 
we obtain the following lower bound on H^"^ : 

H^f^ >~ajdz{Lo+ \Aan, - w]+) ^ x.{^^f + (C^ - 2a) Vc 

-^Jdx |Vx(x)|2 - j^N Isuppxl (^/l^llxlloo + ||Vx||2^ • (3.17) 

The volume of the support of |suppx| = /^Oi enters via the identity JdzUz = 
A^jsuppxl- We further bound [Aariz — t^^]^ < 4an^, and use that 

N N 



J dzUz^XzixiY ^ j dzez{xj)xz{xif' 

i—1 ■ij — 1 

= 2 j dzBz{xj)xz{x,f + N . (3.18) 



l<i<j"<Ar 

Moreover, 



dz9z{xJ)xz{x^f — r (3.19) 
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with 

Z sup (00 * x') (x) . (3.20) 

The final resuh is 

-ff^^^ > ([/ - 2a - %o?Z) Vc - 'ia^N - aNuj 

-f llVxIl^- (^iVlsuppxl (^y^llxlloo + llVxIb^ ■ ^^-2^) 

Note that Z is bounded above by the diameter of the support of x^ which can be 
chosen arbitrarily small. In particular, we can choose the diameter small enough 
such that Sa^Z < U ~ 2a, which leads to a lower bound on Hjj''' that is linear in 
N. This concludes the proof of Theorem [31 □ 

For U — va, v > 2, our lower bound is proportional to a^N for large a. To 
see this, we choose the diameter of the support of x to be of the order 1/a. Hence 
Z ^ a^^, and also cj ~ a by scaling. Moreover, ||Vx||2 ^ a, llxlloo ~ o?^"^ 
and jsuppxl ~ ol~^ , hence the right side of (j3.21l) is of the desired form, namely, 
— const a^N for large a. 

We conjecture that, ioi U — va, v > 2, 

H\f^ > NE'-^'^ (a) ~ C^Na^ for all a > 0, (3.22) 

for some constant depending only on v. For N = 2 this was proved in the 
previous section, but the proof of p.22p for general N > 2 remains an open problem. 



4. The A^-polaron problem: Absence of Binding 



We now return to the question of binding of polarons and prove Theorem [T] 
Because of subadditivity of the energy (|1.8p . E^\a) < NE^^^a) for any N, U 
and a. Hence it remains to prove the reverse inequality. 

We perform a localization similar to that in the two-polaron case, but relative 
to the nearest neighbor. This type of localization is one of the main technical 
ingredients in our proof. As in the bipolaron case, the goal will be to localize each 
particle in a box whose size is of the same order as the distance to the closest 
particles, as long as this distance is not too small. 

Let cpi be given as in (|2.ip - (|2.2p . for some £ > and b > 1. If denotes the 
distance of Xi to the nearest neighbor among the Xj, j ^ i, then 



N 



1= E ni^^--(^^)i' 



(4.1) 



and, by the IMS localization formula. 



(N) 



N N 
-i— 1 j — 1 k 



(4.2) 



We claim that the following bound on the localization error holds. 
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Lemma 5. On the support ofipj.{ti), 

N 



1 52(l-j.) 



3 = 1 k 



{l-llh) 



(4.3) 



with 7 13(7r/2)2 



Proof. Note that ipkitj) depends on Xi in one of two ways. First, through ti when 
j — i, but also through all the tj, j ^ i, where Xi happens to be the nearest neighbor 
of Xj. 

We claim that there can be at most 12 of those Xj. If x is the nearest neighbor 
of both Xj and Xk, then \xj — Xk\ > m.ax{\xj — x\,\xk — x\}, and hence the angle 
between Xj — x and Xk — x is at least tt/B. Think of x as the center of a unit sphere. 
The lines from x to each of these xj 's intersects the unit sphere at certain points 
Pj, whose angular separation is at least tt/S. At each of these points pj we can, 
therefore, construct a unit sphere tangent at pj to the given sphere around x. From 
the packing problem we know there can be at most 12 such spheres. This proves 
the claim. 

On the support of ipj.{ti), 

^2 



E 

k 



ViV?fc(ti) 



< 



(4.4) 



as we have already used in (|2.5p . If Xi is the nearest neighbor of Xj, the same is 
true with ti replaced by tj on the left side, since tj > ti by definition, and the left 
side is easily seen to be decreasing in ti. This concludes the proof. □ 

We now proceed with the one-particle localization as in the two-polaron case, 
localizing particle i in a ball of radius 6^*L centered at Ui, with L < £/(4fe^). More 
precisely, with x given in ()2.7|) . let 



N 



i=l 

where j = (ji, . . • , jw) and u = (ui, . . . , un). We have 

j 



(4.5) 



(4.6) 



and, using Lemma [5l 

. JR3« \ 



N 



i=l 



> 



E 



{N) 



N 
i=l 



76^ 71 



L 1 , . . . , A: 



HVxHi , 

52j-^2 I 



(4.7) 



In analogy with the two-particle problem, the goal here is to show that the integrand 
in this last expression is bounded below by NE'-^\a)\\^jjj^u\\^ which, together with 
()4.6p . implies the conclusion of the theorem. 
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For given j and u, let Bi denote the ball of radius b^^L centered at m. Because 
of our assumption L < £/(4fe^), the balls Bi with ji > 1 do not intersect any of the 
other balls. Let dik denote the distance between ball Bi and ball Bk- 

Recall that the ground state energy can be obtained from the T — )■ oo asymp- 
totics of the functional integral (|1.22p . In the case of relevance here, for states 
having the aforementioned support properties of ipj^-a, the Brownian paths uji in 
the functional integral are confined to the respective balls Bi. In addition to this 
confinement, the paths have the property that at any given time t, the separation 
between any uji(t) and its nearest neighbor among the ujk(t), k i, satisfies the 
conditions according to the support of fj^ . We may relabel the particles such that 
ji = for i < M, and ji > 1 for Af < i < N. The exponential in the functional 
integral is a sum of three pieces, A + B + C, where 

A= > ' > • I 2a / ^^^^^ I , . . ^ - [/ 




\u!i{t) - UJk{t + S)\ Jq \uJi{t) - UJk{t)\ J 

(4.8) 



N 



^ L 2 Jo \LU,it) - UJiit + s)\ ' ^ ' 



k=M+l 



and 



M 



2 -^0 k(i)-c^fc(i + s)| i<,tl'<M-^o \u,{t)-UJk{t)\' 

(4.10) 

We first bound A. For k > M and i k, the distance dik between the balls Bi 
and Bk is nonzero. Since the paths uji and LUk are confined to the balls Bi and Bk, 
respectively, 

N k-1 , 

rri \ ^ \ " I ^Ol U 

- ^ ^ \d~k' d,k + 2L{W^+W^)^ 



k=M+l 4=1 



Similar to 
and hence 



A< -T\U\1 — - 2a 



dik 

k=M+l i=l 

Under the assumption that U{1 — ALb'^/i) > 2a this is negative. We not only want 
it to be negative, however, we also want it to dominate part of the localization 
error, namely. 

Using the fact that miui^^: dik < ^b^'' < Ib"^^^^"^ we can bound 

N N fc-1 

k=M+l k=M+l i=l 
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From (|4.12|) and (|4.13|) . we conclude that 



A + T 



2 \ N 



k=M+l 



as long as 



ALh' 



2£ 

> 2a + — 





7 62 



TT 

2 ' J2 



(4.14) 



(4.15) 



Since we are seeking a lower bound on the energy, i.e., an upper bound on the 
functional integral, (14.15^ then implies that A, together with the localization terms 
coming from M < i < N va (|4.7p . are indeed negative by (14.141) . Therefore these 
terms can be discarded in the functional integral. This concludes the discussion of 
the term A and leaves us with B + C and the remaining localization terms from 
(|4.7|) corresponding to 1 < i < M . 

Since B and C refer to different, now non- interacting, sets of particles (namely, 
M + 1 < fc < and 1 < i < M), we see that the functional integral factorizes. 
The term B is just the exponent in the path integral for {N — M) non-interacting 
polarons, each with its own field. The integral of contributes — M)E'^'^\a) 
to the energy. Henceforth we can forget about B. 

With the aid of our previous linear lower bound for U > 2q!, term C is almost 
as simple as term B. Since U > 2a we can write U — 2a + Z + V with Z and 
V, both positive, to be chosen later. Because of the separation condition for any 
1 < i < M and any time t, the distance between uJi(t) and its nearest neighbor 
among the ujk (t) 's is bounded above by £, and hence 



E 



V 



l<i<k<M 



V 

, . ^ . ^T > —M. 

ij^it) - ujk{t) - e 



The integral of e contributes at least 

AM) 



E)^:i,{a) + M 



V 



(4.16) 



(4.17) 



to the energy. By Theorem |3] this is bounded from below by —MC{2a + Z,a) + 
MV/£, where C{2a + Z, a) is the finite constant implicit in Theorem |31 This term 
is at least 

(the second term being the localization error from (|4.7|) ) provided we take 

y = ^(i?W(a) + ^^^|^ + |^ + C(2a + Z,a)) . (4.18) 
Another way to state this is that U must satisfy 



(4.19) 



for some Z > 0. 

We have thus shown that, for any given j and u, the integrand in the first line of 
(|4.7p is bounded from below by iVi^^^^ (a) ||'(/'j,u|P as long as U satisfies the bounds 
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()4.15|) and (|4.19p (for some Z > 0). In combination with (|4.6p . this concludes the 
proof of Theorem [TJ □ 

There are many parameters in (|4.15p and (|4.19l) : b, £, L and Z. The only con- 
straint on them is b > £ > ALb^ , and each choice gives rise to a computable estimate 
on the critical U . We emphasize that the resulting bound on U is independent of 
N. 

Our bound on the critical value Uc{a) is proportional to a for large a. This 
follows with the choice £ ^ L ^ and b — 0(1), in which case condition (|4.15l) is 
of the form U > const a. With Z ^ a, condition (j4.19p is also of this form for large 
a, since the last term is bounded by const for large a, as shown in the previous 
section. We conjecture that the last term in (|4.19|) is actually bounded by for 
all a, as explained at the end of the previous section, Eq. p.22|) . Assuming the 
validity of p.22[) . our method leads to the bound Uc{ct) < const a for all a > 0. 

5. The Pekar-Tomasevich functional 

5.1. Boltzons for U > 2a.. We shall prove the analogue of Theorem [3] for the 
PT functional. The designation 'boltzons' refers to particles without any symmetry 
restriction. 

Proposition 1. If U > 2a then 



^ ^ ' - \-{0.755)Na^ if U > 3a . ^ ' 

We note that this proves (|1.14p with the constant stated in (|1.16p . 

Proof. We write U = 2{a + 5) with some 5 > Q. Given any V', and hence p^, we will 
use two inequalities to bound the first two terms in the functional (|1.9p in terms of 
p^. The first is the Hoffmann-Ostenhof inequality |T5] (see also Cor. 8.4]), 

YI \V^^iX)\^dX> [ \V^{x)\Ux. (5.2) 



(0.0280) NU^/{U -2a) if2a<U<3a, 



The second is the Lieb-Oxford inequality [19], [201 Thm. 6.1 

3N \Xi — Xj\ 2 



(5.3) 



These two bounds imply that (with ^JfhipjN) 
I-oWl/.w [ /iv7^i2 n a^^rTATl/3^&/A , XAT ff 0(x)2 0(y)2^^ 

(5.4) 



-VY;\^]> f|V0|2-(1.68)C/7Vi/308/3)dx + (57V // f^^^l^^Ldxdy 
N ,/b3 V / ././b3vb3 \x - y\ 



Next, we use Holder's inequality 

dx< ( [ (jy^dx) ( [ (j? dx^ (5.5) 



and Lemma [6] below to conclude that 



(5.6) 
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Using {aP^y^^^abc < ^{aa^ + I3b^ + jc^) for non-negative numbers a,b,c,a, f3,j, 
we see that 

(1.68)C/7Vi/3 J 08/3 da; < J \V^\^dx + SN J J ^^^l^^dxdy 



This, together with ()5.4p leads to the lower bound 



'dx. (5.7) 



While this is true for all U > 2a, the right side is not a monotone increasing 
function of U, which we know the left side to be. Therefore we can say that 
sjj''^ (a) is bounded from below by the maximum value of the right side once U 
exceeds the maximum point, which is U — 3a. This concludes the proof. □ 

Lemma 6. For non-negative functions <f> 

fdx\\—l( ^i^)-^dxdy f \V^{x)\'dx. (5.9) 
/ 47r77R3xR3 \x-y\ Jr3 



Proof. We apply Schwarz's inequality 

{^J^j'dx^ =(b|-V'IW0) 





1 
















p2 







(5.10) 



and recall that \p\ is convolution with (47r|a;|) ^ . □ 

The stability of the PT functional with critical repulsion U = 2a remains an 
open problem. In the fermionic case the answer is affirmative, as was shown by 
Griesemer and M0ller [13]. For the reader's convenience we include the proof here. 

Combining the Lieb-Thirring inequality [22] [20l Cor. 4.1] 

yf \V,ij{X)\^dX>Kq-^/'^ [ p^ixf'^dx. (5.11) 

where K :— |(27r)^/3 — g.lS 5 , with the Lieb-Oxford inequality ()5.3p . one deduces 
via\^]>Kq-^/' [ p4xf/^dx-2ail.68) [ p^{xf'^dx. (5.12) 

The minimization of the expression on the right side under the normalization con- 
straint {x) dx = N leads to the lower bound 



viTm>-^-^e'^Na\ (5.13) 



as claimed in Eq. (|1.17p . 
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5.2. Absence of binding in the Pekar-Tomasevich model. We have just given 
an alternative proof of the fact that there is stability of matter in the PT model 
when U > 2a. Now we discuss the other part of Corollary [1] that is, the absence 
of binding for large U. 

The first step is to linearize the problem. The variational problem for the PT 
functional can, equivalently, be written 



4^)(a) = mf ^(^ 



U n^ul A: I \i^?dX = l\ (5.14) 



where the iV-particle Hamiltonian H^^l is defined to be 

+ „/■/■ ^Mf^Li^i,, (5.15, 

JjR^xmi \x - y\ 

We proceed as before, by localizing particles into individual boxes, with sizes de- 
pending on the distance to the nearest neighbor. In each localization region we 
obtain a lower bound on the energy of a given -0 by minimizing over cr, which yields 
the PT functional for the localized ip. (In other words, we linearize, localize, and 
de-linearize. If we had not followed this route and tried to deal with the quartic 
term directly, the resulting expressions would be much more complicated.) 

Consider first the case N = 2. In the region j = 0, we need a lower bound on 
£q (a). We could use Lemma [2] above, but it is simpler, and indeed more accurate, 
to use 

f^''(a) = 2£:('^(2a) =8£^''\a) = 2f^'^(a) - 6 • (0.109) . (5.16) 

The first equality follows from the linearization ()5.14|) since the ground state of 
()5.15p for U — is SL product function for every cr; the second equality follows by 
scaling. 

For the regions j > 1, we obtain the PT functional for 2 particles localized in 
disjoint balls. The proof of the corresponding lower bound to the energy, analogous 
to Lemma [TJ is obvious, bounding the attractive energy using the smallest possible 
distance of the particles. Alternatively, one can use Lemma [H For j > 1 the 
resulting condition on U is thus the same as in the proof of Theorem [5] Our 
improved estimate in the j = region leads to the bound Vc < 29.4 (to be compared 
with the bound C < 2 • (26.6) = 53.2 for the Frohhch polaron). 

We can similarly analyze the A^-particle problem. For the particles with j > 1 
the bounds are exactly the same as before, except that functional integrals are not 
needed in the derivation. For the particles with j = the improved stability bound 
(|5.ip (or (|1.17p for fermions) is used, and hence the final condition for the absence 
of binding will be a lower value of U than that for the Frohlich Hamiltonian. □ 



Appendix A. Proof of Lemma [3] 



We assume that the confining balls Bi , B2 are each of radius i?, and that the 
balls are of distance d = inf{|a;i — X2\ ■ xi £ Bi, X2 £ B2} from each other. In the 
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following, let 

be normalized annihilation and creation operators, with x a point in configuration 

(2) 

space. In terms of these operators, the particle- field interaction term in is 
given by 



Va(j>{x) = \ dy— (A.2) 

TT'^/^ J \x — y\^ 

and the phonon field energy by 

Hf = J dxa''{x)a{x). (A.3) 

Fix a plane midway between the two balls and perpendicular to the line between 
their centers, and let Si and 5*2 be the resulting half-spaces with Bi C Si and 
B2 C 5*2. Then we have the identity 

(2) _ 2 /■ aHy) + aiy) 

^u=o ^ Pi ^ s/? / y I 12 

+ 



dy(dHy)- 3/,,^ 12 ) (^(y) 

2 ^ f ^ a''{y)+d{y) 
P2 Tn / "■y~i 12~ 



a f dy ^ f '^y 



^ 'si |a;2 - y\^ 7r3 |a;i - y\^ 



(A.4) 



Define 



00.2 (2/) = — hyeS'i (A.5) 

V TT'^'-^\x2 - yr J 

and analogous expressions for dl.^{y), and for aj;^ (y), aj.^ (y) with y £ S2- In terms 
of these operators, the identity (jA.4p becomes 

H'Jlo-Pl-^f dy^ll^+f dyaUy)d.M 



71-3/2 1^;^ _ y|2 

2" (A.6) 



7r2|a;i-7T,| 7r2|x2 • n| l^i— X2I' 

where \xi ■ n\ is the distance between Xi and the dividing plane, i = 1, 2. (The inte- 
grals in the last line of Eq. (|A.4p can be done explicitly via cylindrical coordinates, 
resulting in the first and second terms in the last line of this identity; the last term 
is an integral of \xi — y\~''^\x2 — yl"^ over all of and is readily computed to be 
the Coulomb attraction term.) 

We can give a lower bound on expectations of the right side of the first line of 
this last equation (|A.6p . assuming the first particle is indeed confined in the ball 
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Bi. Let Ksj^^x2 be the one-particle operator of the first fine, 



dy 



\xi - 



(A.7) 



which we regard as acting in the Hilbert space (E) -^Sn the latter factor 

being the Fock space associated with the phonon variables y G Si; the operator 
is a function of X2 G S'2 • Note that pi commutes with 0,^2 {y) and its adjoint and 
that (y) and aj,^ {y) satisfy the canonical commutation relations. The operator 
function Ks2.xi is defined analogously. 

Fix X2, let ip he a state in L'^{M.^)(E)J-Si supported in xi G and then consider 
a product state "if = ip^^ G L^(R^) (g) T where $ is a coherent state of the phonon 
variables corresponding to y G 52 such that 



t3/2| 



y\ 



\^), ye S2. 



(A.8) 



Here, we take Xc S -S*! to be on the line passing through the centers of the two balls 
and of distance d/2 + 2R from the dividing plane (i.e. as remote from the dividing 
plane as possible but on the surface of Si). For such a state ^P, we have that 

= {'4'\Ksi,x2\ip) - ^ / dy(tp — — ^ 



\Xl 



X2m — T 



inf 



-yl' 

dy 



a 



dy 



S2 



\Xl 



y\ 



7r3(d/2 + 2i?) 



(A.9) 



The integral in the infimum is seen to have no critical points for xi in the interior 
of Bi and so attains its minimum for xi on the boundary of Bi. The integral 
can again be written using cylindrical coordinates and the angular integration per- 
formed explicitly. One then writes the integrand of the resulting double integral 
just as a function of zi, say, where xi — {ri,9i,zi) in cylindrical coordinates, and 



where (zi -I- d/2 + RY 



R . Minimization of the integrand in this double 



integral regarded as a function of zi is tedious but straightforward, the minimum 



The integral /^^ dy \ 



occurring at xi = 
computed above for Eq. (IA.4p . Thus, we obtain 

E'^^Xa) < {iJi\Ks„x2\i^i) 



y\ 



is equal to 7r/(c?/2 + 2R) as 



a 



(A.IO) 



TT^{d/2 + 2R)' 

Of course the second line of Eq. (jA.6p is handled similarly. 

By this last inequality (jA.lOp and Eq. (jA.6[) . we have that for any state "if with 
electron support in Bi x B2, 



1|*> 
1 



{ms2,.. 
1 



11*) 
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Noting that \xi ■ n\ and \x2 ■ n\ are at least d/2, we have that 

which is the claim of Lemma |3l □ 
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